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ABSTRACT 



Context. Though only a handful of extrasolar planets have been discovered via direct imaging, each of these discoveries 
had tremendous impact on our understanding of planetary formation, stellar formation and cool atmosphere physics. 
Aims. Since many of these newly imaged giant planets orbit massive A or even B stars we investigated whether giant 
planets could be found orbiting low-mass stars at large separations. 

Methods. We have been conducting an adaptive optic imaging survey to search for planetary- mass companions of young 
M dwarfs of the solar neigbourhood, to probe different initial conditions of planetary formation. 

Results. We report here the direct imaging discovery of 2MASS J01033563-5515561ABb, a 12-14 Mj up companion at a 
projected separation of 84 AU from a pair of young late M stars, with which it shares proper motion. We also detected 
a Keplerian-compatible orbital motion. 

Conclusions. This young L-type object at planet/brown dwarf mass boundary is the first ever imaged around a binary 
system at a separation compatible with formation in a disc. 
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Stamatellos et al.||201l[ ) 



The discovery of hundreds of extrasolar planets in the last 
20 years has radically modified our understanding of plan- 
etary formation. Though radial velocity and transit detec- 
tion methods have proven by far the most prolific, the few 
planetary-mass companions which have been discovered by 
direct imaging have provided very challenging constraints 
for formations model s, especially the core-accretion model 
(Pollack et al. 1996| ) that is prefered to explain the for- 
mation of Sola r Syste m planets. 2M1207B, discovered by 
( |2004[ ) , with a mass-ratio of 20-25% is too 
respect to its primary to have formed by 



Circumbinary planets, such as Kepler-16 ABb (Doyle et al 



Chauvin et al. 
with 



massive 



core accretion, while most of HR8799 (Marois et al. 2008) 



2011) are even rarer and provide peculiar constraints on 
planetary formation scenarios, notably on the influence of 
binarity on planet-forming discs. 

We present here the discovery of 2MASS J01033563- 
551556lABb, hereafter 2MASS0103(AB)b, a unique a 12- 
14 Mj up substellar companion to a late M dwarf binary 
system. 

2. A 12-14Mj up companion orbiting around a 
young late M binary system 



would be very di fficult to form in situ by core-accretion. 2.1. Observations and data reduction 



Only /3-Pictoris b ( Lagrang e et al.|2010 ) fits relatively well 
with the core-accretion s cenario. Also, several imaged sub 
stellar companions 

"1, 



et al. 
and 



2008; (Carson et 
ibated- 



C hauvin et al. 2005 Lafreniere 
2012[ ) straddle the arbitrary - 
13 Mj up planet/ brown dwarf boundary. For 
most of these massive planets (or light brown dwarfs) the 
formation mechanism, stellar or planetar, is s t ill debated 
( Luhman et aT1|2006| |Bate||2009[ |Rafikov|[MT| |Boss||2011 



Send offprint requests to: P. Delorme, e-mail: 
Philippe.Delorme@obs.ujf-grenoble.fr. Based on ob- 
servations obtained with NACO on VLT UT-4 at ESO-Paranal 
(runs 090.C-0698(A) and 70.D-0444(A). 



We imaged 2M0103 in November 2012 (run 090.C- 
0698(A)), in L' band as one target of our NACO s urvey for 
planetary c ompanions to young nearby M dwarfs ( Delorme 
et al.|2012[ ). We used NACO infrared wave-front sensor and 
observed in pupil tracking (only 12° of rotation) and cube 
mode in V , and our follow-up observations in JHKs on the 
same night used field tracking. Table [T] shows the details of 
our observations. 

The target star was resolved as a low contrast, 0.25" 
binary on these raw images and an additional source was 
identified at -1.8" at the north west of 2M0103A. In order 
to measure the proper motion of this source, we retrieved 
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Table 1. Summary of the NACO (VLT-UT4) observations 
of 2M0103AB (RA=01:03:35.63; Dec=-55:15:56.1). 



UT Date 


Filter 


Exp. time 


Comments 




U 


32x200x0.2=1280s 


Seeing: 


2012-11-25 


K s 


8x20=160s 


0.7"-0.8" 




H 


8x20=160s 


Airmass: 




J 


4x5=20s 


1.16-1.25 


2002-10-28 


H 


5xl0x2=100s 


Archive data 



J 



O 



N 



0.3" 



Fig. 1. 2MASS0103(AB)b in November 2012, with NACO 
in L' band. The green arrow shows the position of the com- 
panion in 2002. The light-blue circle identifies the expected 
position of the companion if it had been a background 
source. 



ESO archive NACO #-band images of 2M0103, obtained 
in October 2002 (run 70.D-0444(A)). These early images 
were acquired in field tracking and with poor adaptive 
optics correction. We stacked the best 50% of the frames, 
for which the central binary was resolved, totalling 100s 
exposure time on target. 

We used t he IPAG-ADI pipeline as described in |Delorme| 



et al. (2012) to reduce the frames (bad pixel interpola- 
tion, flat, recentring, derotation and stacking). Although 
both the secondary component and the companion appear 
clearly after a si mple stack of all ex posures (see Fig{T]), we 
performed ADI (Marois et al. 2006) and LOCI (Lafreniere 



et al.|[2007 ) star subtraction procedures to detect eventual 



other companions. None was detected was detected down 
to a contrast of ~7.5 magnitudes at 0.5", resulting in a de- 
tection limit of ~2.5 Mj up at 25 AU for an age of 30 Myr 
(see discussion below). 



Tucana-Horologium (TH A), Columba, Carin a, Argus and 
AB Doradus (ABDMG) ( |Torres et atl[2008| . The details 
of this analysis will be presented in Gagn et al. (in 
prep.), but the principle is to identify promising candidate 
members to these moving groups using astrometry, proper 
motion and photometry from a correlation of 2MASS and 
WISE catalogs, with a modified version of the Bayesian 
analysis described in Malo et al. ( |2013[ ). One of the first 
robust candidates identified in this project is 2MASS 
J01033563-5515561, which we have followed with GMOS-S 
at Gemini South to obtain the optical spectra. This 
spectrum matches a M5.5/M6 spectral-type and shows 
strong H-alpha emission at 656 nm, with an equivalent 
width of 10.23±0.55A No nearby X-ray so urce was found 
in the ROSAT archive ( |Voges et al. 1999), indicating the 
target is not a strong X-ray emitter. In parallel to this, 
we have obtained a trigonometric distance, of 47.2±3.1 pc 
for this object (A. Riedel, private communication, using 
the CTIO 0.9m through the CTIOPI program, using 49 
i?-band images taken on 11 nights between October 26th, 
2007 and November 13 t h, 2012, and reduced using methods 
from [Jao et al.| ( |2005| ); [Riedel et al.| ( |2011[ ). The complete 
parallax analysis for 2M0103, together with many other 
objects, will be published in Riedel et al., in preparation. 
During the NACO runs described earlier, we have also 
noticed the primary is in fact a binary with a flux ratio of 
0.8 in the L' band. Taking into account this binarity and 
the trigonometric distance, we find Bayesian probabilities 
of 99.6% and 0.4% for membership to THA and ABDMG 
respectively. The field hypothesis has a probability of 
10 -14 . 2M0103AB is therefore a strong candidate member 
of the Tucana-Ho rologium association, aged ~30 Myr 
(|Torres et aLpX)8 ). 

We must stress that those probabilities are not absolute 
ones in the sense that even a sample of candidates with a 
100% Bayesian probability will contain a certain number of 
fa lse-positives . Follo w-up observations of robust candidates 
in Malo et al. ( 2013[ ) have shown that the false-positive rate 
is 10% for candidates without parallax in THA. Though 
the membership analysis in our study is not exactly the 
same, the risk of a false positive is very low, especially 
because we do have a parallax measurement, meaning that 
2M0103AB is very probably a bona-fide member of THA 
We will assume in the following that the 2M0103 system is 
aged 30 Myr. 



2012 



Ac cording to BT-Settl 2012 isochrones (Allard et al 

a distance of 



Baraffe et al. 



2003| ), and assuming 
47.2±3.1 pc and an' age of 30 Myr, 2M0103AB is a low 
mass binary with masses of [0.19;0.17]±0.02 M Q for [A;B] 
respectively, see table [2] The projected separation between 
A and B was 0.26±0.0F in 2002 and 0.249±0.003" in 2012. 
The projected distance was around 12 AU at both epochs, 
but the position angle changed significantly, from 71.2° in 
2002 to 61.0° in 2012 



2.3. Proper motion analysis: a bound companion 



2.2. Host star properties 

The primary star 2MASS J01033563-5515561 was identi- 
fied as part of a survey designed to identify new, later than 
M5 candidate stars and brown dwarfs to the young, nearby 
moving groups and associations Beta Pictoris, TW Hydrae, 



During our November 25th, 2012, L' band NACO observa- 
tions of 2M0103 (run 090.C-0698(A)), we identified a can- 
didate companion with a separation of 1.78" and a position 
angle of 339.3° from the primary 2M0103A. Even if contam- 
ination by background objects is rel atively low in L' ban d 
compared to shorter wavelength (see Delorme et al.|[2012 ), 
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Table 2. Host system absolute magnitudes compared 
with BT-Settl isochrones at 30Myr absolute magnitudes 
(2MASS for JHK and NACO for L'), 



Filter 



Mj 



M 



H 



M. 



Mr 



2M0103A 
2M0103B 
Model O.2M 
- 0.175M Q 



7.36±0.05 
7.56±0.05 

7.31 

7.51 



6.78±0.05 
6.98±0.05 

6.75 

6.95 



6.44±0.05 6.04* 

6.64±0.05 6.24* 

6.50 6.1 

6.70 6.3 



Notes. * Since no calibrated photometry is available in L' , these 
magnitudes are derived from the modelled Ks — L'. 

Table 3. Separation (Sep.) and position angle (PA) of the 
companion, with respect to 2M0103A and to the center of 
mass of the binary. 





2002-10-28 2012-11-25 


Sep. from 2M0103A(") 
PA from 2M0103A(°) 


1.682±0.015 1.784±0.003 
341.7±0.05* 339.8±0.01* 


Sep. from baryccnter(") 
PA from barycenter(") 


1.718±0.015 1.770±0.003 
338.0±0.05* 336.1±0.01* 



Notes. * The error in position angle refer to the relative error 
between both epochs. The absolute error, dominated by system- 
atic uncertainties in the position of the reference stars in theta 
Ori., is ±0.4°. 



the probability that this companion was a contaminant was 
high. However, the companion was redder in Ks — L' than 
its late M host system, meaning that the companion is of 
even later spectral type, considerably decreasing the like- 
lyhood of the contaminant hypothesis, but not incontro- 
vertibly proving that the companion is bound. A definite 
proof of companionship was however provided thanks to 
archive images taken with NACO on October 28th, 2002. 
Using this 10 years time base, we could determine that 
2MASS0103(AB)b unambiguously (contamination proba- 
bility <0.001%, taking into account parallax motion) shares 
the proper motion of 2M0103AB (see Figfl]) and is there- 
fore a bound companion. The relative astrometry at each 
epoch is shown on table [3j 

2.4. Companion properties 

Our current information about the physical properties of 
the system relies on the J, H, K and L' photometry and 
astrometry from our November 2012 NACO run, as well 
as October 2002 NACO archive images. The resulting 
Moffat-fitting photometry is J=15.4±0.3, iJ=14.2±0.2, 
if s =13.6±0.2 and L'= 12.6±0.1. We emphasize that all 
these measurements were derived relative to 2M0103A, 
which introduces significant systematic errors (hence 
the large errors bars). The absolute magnitudes of the 
companion, for a distance of 47.2 pc, are shown in table [4] 
and are compatible with 2MASS0103(AB)b being a 12 to 
14 Mj up companion orbiting at 84 AU around the young 
low-mass binary 2M0103AB. Note that 2MASS0103(AB)b 
very red colours in JHKs do not match the colours of field 
objects of similar absolute magnitudes and are similar to 
known, early L-type, young planetary-mass objects, thus 
independantly confirming the youth of the system (see also 
Fig.§. 



It is to be noted that the age of THA is not perfectly known 
and the dispersion of the age estimations of individual 
stars in THA span the 20-50 Myr range (Zuckerman & 
Webb||2000| |Torres et al.||2000[ ). If we assume an age of 20 
Myr, 2MASK0l03(AB)b would be a 12-13 M Jup planet, 
while it would be a 14-15 Mj up brown dwarf if we assume 
an age of 50 Myr. In spite of the naming change, the 
physical differences in mass estimates for the 20 Myr 
and the 5 Myr hypothesis ar e much smaller than those 
derived in Marois et al. (2010) for HR8799bcde planets in 
the same age range. An explanation is that objects more 
massive than ~10 Mj up undergo some deuterium burning 
in this age range, somewhat compensating cooling down 
mechanisms. However, since there is currently no robust 
independent mass constraint for any i maged exoplanet 
(save d to some extent for /3Pic. b, see Lagrange et al. 
2012), it is probable that the systematic uncertainties 



coming from substellar models inaccuracies are larger 
than those arising from age uncertainties. For the sake 
of comparison with other substellar companions found in 
associations of the same mean age of 30 Myr, we assume 
in the following that 2MASS0103(AB)b is a 12-14 M Jup 
object aged 30 Myr. 

A possible analog, if confirmed as bound, would be the 
substellar object located at a projected sep aration of 



1100 AU from the binary system SR12AB (iKuzuhara 
eFaI7|[2qTT] > . The properties of 2MASS0103(AB)b (mass of 
12-14 M Jup , age of 30 Myr, colours, projected separation 
of >50 AU) a nd observed colours are also much like 
AB pic b (K1V Chauvin et al. 2005) or k Andromeda b 
( Carson et al.|201^[ ). The properties of the host systems are 



however quite different. While k Andromeda is a massive 
B star (~2.5M , mass ratio of -0.5%), 2M0103AB is a 
close binary system composed of 2 late M dwarfs, whose 
combined mass is — 0.36M Q , resulting in a mass ratio of 
approximatively 3.6% for the system. 



The position of the companion at each epoch was de- 
rived by Moffat-fitting and the orientation of the detector 
was calibrated using NACO calibration images of theta Ori. 
obtained close in time of the science images. As shown on 
table[3j the relative astrometry is accurate enough to detect 
the orbital motion of the companion around the center of 
mass of the system, with a projected motion of 77±15 mas 
over ten years. The corresponding velocity at 47.2 pc is 
1.7±0.3 km.s -1 . The Keplerian velocity, assuming a circular 
orbit of 84 AU around the 0.36M Q system is 1.96 km.s -1 , 
corresponding to a period of 1280 years and is fully com- 
patible with our measurement. It is to be noted that the 
secondary and the companion rotate in the same direction 
and that their observed orbital motion can be compatible 
with a face-on orbit but not with an edge-on one. 



3. A challenge for stellar and planetary formation 
theories 

2MASS0103(AB)b has a companion mass to host system 
mass ratio of ~0.036, which is too low to match known 
low mass multiple systems (See Fig{3]and also Allers et al., 
2007), but still higher than most star-planet systems con- 
firmed so far. Systems with similar mass ratio, but almost 
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Fig. 2. J — Ks versus K$ — L colour-colour diagram show- 
ing 2MASS0103(AB)b together with other known planetary 
and brown dwarf companions to youn g stars (taken from 
Bonnefoy et al.|2013||Bailey et al.l2013[ ). The symbols with- 
out error ba rs show the colour of he ld M, L and T dwarfs 
(taken from Golimowski et aL||2004 ). 



Table 4. Companion absolute magnitudes compared with 
BT-Settl isochrones at 30Myr and 5Gyr (field hypothesis) 
absolute magnitudes, and other known companions at the 
planet/brown dwarf mass boundary. 



one order of magnitude more massive, with 55-95 Mj up 
brown dwarfs or very low-mass star s orbiting massive ( 1-2- 
2M©) stars, have been identified by Janson et al. (2012 1. In 
case formation mechanisms would simply scale with mass, 
the same processes could be at work to explain the forma- 
tion of these very different kind of objects. If we keep to 
more similar low-mass star systems shown on Fig. [3j the 
mass ratio of 2MASS0103(AB)b is very close to those of 
DH tau B (8-22 M Jup , separation of 330 AU) and CHXR 



73b (7-20 M J? , sep aration of 210 AU) ( |Itoh et al.||2005 
Luhman et al.||2006[), but its projected separation is much 
smaller. Luhman et al. ( 2006 ) state that neither DH tau B 
nor CHXR 73 b could be formed in situ by core-accretion 
or disc instability because of the very large separation from 
their host stars, and the same holds for the 1100 AU can- 



0.35 



0.30 



0.25 



o 



, 0.20 



0.15 



0.10 



0.05 



0.00 



.Sr. 



• 


Imaging 


• 


Transits 


• 


Radial velocity 


• 


Micro-lensing 


• 


DH Tau b 




CHXR 73b 


• 


2MASS0103ABb 



o.o 



0.2 



0.4 0.6 
star mass / solMass 



0.8 



1.0 



Fig. 3. Substellar companion to stellar host mass ratio 
versus stellar host mass, sorted by their discovery tech nic, 
from the exoplanet.eu database (Schneider et al. 2011) 
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didate companion to SR12AB ( Kuzuhara et al.[2011 1. The 
case is different for 2MASS0103(AB)b, at a separation of 
only 84 AU. At such separations, a formation in a gravi- 
tationally instable primordial circumbinary disk would be 
fully compatible with plane tary formati on by gravitational 
instabilities, as de scribed by Boss (2011 ). Howev er, this sce- 
nario is discussed: Dodson-Robinson et al.| (2009 1 claim that 
objects formed by disc instabilities around M-dwarfs should 
have ~10% of the mass of the host system meaning that 
2MASS0103(AB)b would not b e massive enough for such 



a scenario, while other studies ( Rafikov 2009 Stamatellos 
et al.|[2011 ) find that such low-mass discs cannot fragment 
at all. Simultaneous formation and ejection of the 3 com- 
ponents in the massive disc of a more massive orginal host 
st ar is plausible, in a scenari o akin to what is described 

but 



Stamatellos & 



Whitworth 
with 



^ ([2009}, but the central bi 
masses of 0.17 and 0.19 M Q 



arc 



nary components, 

more massive tha n most objects formed in [Stamatel los fc| 



Whitworthj ([2009J simulations. 

A planetary formation scenario by core-accretion (e.g. 



[ Kennedy fc Kenyon|[2008[ |Mordasini eraL[|200"9] |Rafikov| 
|201ip can very probably be excluded for several reasons. 
First, the separation is too large for a formation in situ. 
Second, the companion has ^3.6% of the mass of its host 
system, which is of the order of magnitude of the maxi- 
mum total mass of the protoplanetary disc from which core- 
accretion planets are formed. Finally, such a 12-14 Mj up 
companion would be a very rare occurence, according 
to the core-accr e tion p lanetary mass function derived by 



Mordasini et al. (2012) 



A purely stellar for mation mode by turbulent core 
fragmentation (see e.g. [Padoan fc 
20091 |Hennebelle & Chabrier||2011|) 
Tb 



Hennebelle & Cha 
tETcase 2MASS0103(AB)b would" 



Nordlund 


2002 


Bate 


} is plausi 


Die, and in 



c an extreme case of 
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hierarchical triple stellar with a third component in the 12- 
14 Mj up mass range. However, a stellar formation scenario 
would necessitate that cores can naturally fragment into 
such low mass objects, without requiring any ejec tion from 
the the accretion reservoir (such as des cribed in Reipurth 



fc Clarke 
be difficu 



2001 Bate & Bonncll 2005), because it would 



t to starve the accretion of the third component 



without also stopping accretion on the central binary. From 
hydrodynamica l simu l ations of stellar formation by cloud 
fragmentation, 



Bate 



(2012) claims that "brown dwarfs 
with masses <15 Mj up should be very rare", implying that 
formation by direct core fragmentation of a 12-14 Mj up 
object such as 2MASS0103(AB)b would be possible but 
uncommon. 

In any case, the discovery of 2MASS0103(AB)b brings 
most current stellar and planetary formation theories to 
their limits while others, such as core-accretion, can prob- 
ably be excluded. The very existence of such a peculiar 
system therefore provides a very valuable test case against 
which current and future stellar and planetary formation 
theoretical models can be tested. 
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